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Abstract. We investigate spin-dependent decay and intersystem crossing in the 
optical cycle of single negatively-charged nitrogen-vacancy (NV) centres in diamond. 
We use spin control and pulsed optical excitation to extract both the spin-resolved 
lifetimes of the excited states and the degree of optically-induced spin polarization. 
By optically exciting the centre with a series of picosecond pulses, we determine the 
spin-flip probabilities per optical cycle, as well as the spin-dependent probability for 
intersystem crossing. This information, together with the indepedently measured decay 
rate of singlet population provides a full description of spin dynamics in the optical 
cycle of NV centres. The temperature dependence of the singlet population decay rate 
provides information on the number of singlet states involved in the optical cycle. 
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1. Introduction 

Nitrogen-vacancy (NV) centres in diamond are well-defined quantum systems in the 
solid state, with excellent spin coherence properties [TJ. Even at ambient conditions, 
NV centres have successfully been used in the field of quantum information processing 
El m El El El E], magnetic sensing P QUI [EEJ H21 03] and photonic devices [Ml E51 
[TBI IT71 HU [19] . However, despite the rapid experimental progress the understanding of 
the optically-induced spin-dynamics of the NV centre is still incomplete. In particular, 
the spin-dependent intersystem crossing (ISC) rates as well as the number of singlet 
states involved in the optical cycle are still debated. These parameters are responsible 
for optical spin initialization and readout, and are important for a correct estimation 
of photon emission rates. We extract these values by a series of room-temperature 
experiments, where we perform spin-resolved fluorescence lifetime measurements using 
picosecond optical excitation pulses. The lifetime of the singlet manifold is measured 
by analyzing the initial fluorescence rate for consecutive microsecond optical pulses 
with variable delay. The temperature dependence of this lifetime yields insight into the 
number of involved singlet states. 

2. Experimental setting 

We investigate individual NV centres contained in a high-temperature high-pressure 
(HTHP) grown type Ha diamond sample from Element Six ((111) oriented). The 
sample is studied in a scanning confocal microscope setup operated at T = 10 . . . 300 K. 
Spin control is achieved via microwave (MW) fields applied to an Au-waveguide that is 
lithographically defined on the diamond surface [6J, E] • For optical excitation we use a 
continuous-wave (CW) laser at A = 532 nm, equipped with an acousto-optic modulator 
(AOM) with 20 ns rise-time, as well as a frequency-doubled diode laser at A = 532 nm 
with a pulse length of 62 ps (max. pulse energy: 25 nJ) and variable repetition rate. For 
photon detection we use an avalanche photo diode in the single-photon counting regime 
with a timing jitter of 450 ps. Time-resolved data are acquired using a time-correlated 
single photon counting module with a jitter of 10 ps, using a bin size of 512 ps. An 
arbitrary waveform generator (channel-to-channel jitter < 100 ps) is used as timing 
source of the experiment. 

3. Model 

The photo-dynamics of the NV centre [figure [T]( a)] are determined by six electrons, which 
in the ground state form a triplet occupying an orbital of 3 A 2 symmetry. The centre can 
be excited via a dipole-allowed transition to a 3 E triplet state. This level also has a spin- 
dependent probability to undergo ISC to a series of singlet states [20J. We use a five- level 
model to describe the spin dynamics of the NV centre [figure [T^b)]. Spin-spin interaction 
splits the 3 A 2 ground state by D GS = 2.87 GHz into a state with spin projection m s = 




Figure 1. (a) Lattice structure of the NV centre: a substitutional nitrogen atom N 
next to a vacancy V in the diamond lattice C. (b) Level structure of the NV centre: we 
consider spin-conserving (fei, &42) and spin-flip (/C32, k&i) transitions between triplets 
(states with spin projection m a = ±1 are merged into states |2), |4)). Spin-dependent 
ISC rates connect triplets to the singlet states (summarized as |5), and described by 
rates £35, &45, &51, £52)- 

(|1)) and a doublet with m s = ±1 (summarized in |2)). Correspondingly, the excited 
3 E state is labelled |3) {m s = 0, associated with a lifetime T 13 ) and |4) [m s = ±1, 
lifetime T 14 ), and split by D ES = 1.43 GHz [211 [22]. Two singlet states with a splitting 
of AE = 1.189 eV have been identified experimentally [23[ El], but recent theoretical 
studies [25] and also data obtained in this work suggest the presence of a third singlet 
state between 3 A2 and 3 E. For the analysis of spin dynamics, however, we summarize 
the singlet states in |5), and the corresponding lifetimes are summed and denoted as 
Ti i5 . Rates from state \m) to state \n) are denoted by k mn , and we only consider 
relaxation rates indicated in figure H](b). Population in state \n) is denoted by P n , and 
the spin polarization in ground and excited state are denoted by Pes = Pi/ (Pi + P2) 
and P E s = P 3 /(P 3 + P 4 ). 

4. Spin-dependent lifetime 

Pulsed optical excitation and time-resolved detection of fluorescence provides a simple 
and direct way to determine the lifetime of the excited state in an optical transition. 
If the excitation pulse is short compared to the lifetime 7\, the detected time- 
resolved fluorescence (averaged over many excitation cycles) decays exponentially / oc 
exp(— t/Ti). However, if the system under consideration is excited into a mixture of n 
excited states with different lifetimes, the detected fluorescence decays according to a 
multi-exponential function / oc J2 n a n exp(— t/T\ %n ). 

This situation is present in the case of NV centres in diamond, where the excited 
state is composed of a spin triplet. Population in this state decays radiatively to the 
triplet ground state with identical oscillator strength for the different spin projections 
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Figure 2. (a) We first polarize the spin by applying al.3/is laser pulse at A = 532 nm. 
After 800 ns we turn on a MW field at 2.87 GHz for a variable duration. 200 ns after 
begin of the MW we excite the NV center by a 62 ps laser pulse at A = 532 nm and 
measure the time- resolved emission. The first 300 ns of the subsequent polarization 
pulse is used for spin readout. The experiment was performed at T = 300 K. 
(b) Fluorescence decay as function of MW burst duration (fluorescence counts are 
encoded in a logarithmic colour scale). The decay time oscillates with the m s = 
amplitude, as confirmed by (c) conventional spin readout, (d) Fluorescence decay 
curve, integrated over all applied MW burst durations. A fit using a bi-exponential 
function, yields the lifetime of states |3) and |4). A single-exponential function (shown 
for comparison) cannot accurately fit the experimental data, (e) Degree of spin 
polarization Pes = P3AP3 + Pi) with no MW applied and after a 20 ns MW pulse, 
obtained from the relative amplitudes of a bi-exponential fit. 



[26J. Because of spin-dependent ISC [20], the effective lifetime of the excited state 
is significantly different for states |3) and |4) . In literature [221 EEJ [27] the NV 
centre's lifetime generally is obtained by fitting the time-resolved fluorescence to a 
single-exponential decay, leading to values of Ti i3 ps 12 — 13 ns and Ti j4 ps 8 ns in 
bulk diamond. Since optically induced spin polarization in NV centres is limited 
[28| [29] . in fact we expect such a lifetime measurement to yield a bi-exponential 
decay curve with time constants set by the sums of rates out of states |3) and |4) 
[T 1)3 = 1/(^31 + /C32 + ^35) and T 1;4 = l/{k±\ + k i2 + £45)], and an amplitude ratio set 
by the initial spin polarization. Such a bi-exponential decay has been observed in [29J, 
where a polarization of Pes — 0.84 ± 0.08 % has been obtained. These data are based 
on fluorescence lifetime measurements with MW spin manipulation in the excited state, 
where the duration of the MW pulse was neglected, and the pulse was assumed to be 
perfect. 

Here we present a simple and reliable way to obtain the spin-dependent lifetimes, 
without any assumptions on the MW pulse and spin polarization. For that purpose 
we drive Rabi oscillations in a conventional fashion, i.e. we apply a 1.3 fis long off- 
resonant laser pulse to polarize the electron spin, followed by a MW pulse of variable 
duration, resonant with the zero-field splitting of D = 2.87 GHz. After the MW, we 
apply a ps laser pulse. The fluorescence as function of MW pulse duration is shown 
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in figure |2](b). For comparison we also plot the result of a conventional spin readout 
[figure [21(c)], i.e. the intensity of the first 300 ns of the subsequent polarization laser 
pulse [30]. The data in figure 121(b) clearly reveal the oscillations in decay time which 
are in phase with oscillations of the electron spin. To accurately fit these data to a bi- 
exponential decay we sum over all MW pulse durations, which gives similar contribution 
from each spin state. From the fit we obtain the two time constants T 13 = 13.7 ± 0.1 ns 
and T 14 = 7.3 ± 0.1 ns [figure 121(d)]. By using these values as constants for a fit to the 
individual decay curves, we can determine the relative contributions of states |3) and 
1 4) to the minima and maxima of the Rabi oscillation data from the relative amplitudes 
of the two exponentials. For this particular centre, we find PES,max = 72.1 ± 0.9% 
and PEs,min = 12.2 ± 0.5% [figure [21(e)]. The value of Pcs,max m ay be larger due to a 
nonperfect spin conservation in optical excitation (see section [6]). 

We note that for a spin S = 1 system with degenerate levels m s = ±1 (i.e. for 
NV centres in absence of magnetic field), the effect of resonant MW driving starting 
from a pure \m s = 0) state is to cause coherent oscillations between \m s = 0) and the 
symmetric superposition \m s = +1) + \m s = —1). In a more realistic scenario we need 
to consider an only partially polarized state, where the density matrix subspace spanned 
by \m s = +1) and \m s = — 1) has equal population of the symmetric and antisymmetric 
superposition state. Only the symmetric state will be transferred back to \m s = 0), so 
the effect of a half-oscillation is to transfer the full m s = population into m s = ±1, 
while only half of the incoherent population in m s = — 1 and m s = +1 is transferred 
back into m s = 0. As a consequence, when driving Rabi oscillations, the maximum 
population in m s = will only reach half the value of the maximum population in 
m s = ±1 [figure [2](e)]. 

In summary, the presented method allows for accurate determination of the lifetime 
of the pure states |3) and |4), without the need of assumptions on the quality of spin 
manipulation. The knowledge of these lifetimes can then be used to quantify the spin 
polarization. 

5. Temperature dependence of singlet decay 

An important parameter for the spin dynamics of the NV centre under optical excitation 
is the time population spends in the singlet manifold before it decays back to the triplet 
ground state. This time scale is long compared to the excited state lifetime, and since 
decay into the singlet states (ISC) is favoured for m s = ±1, this leads to reduced 
fluorescence for m s = ±1, a fact which is routinely used for non- resonant spin readout 

The temperature dependence of this decay rate yields information about the energy 
splitting involved in this relaxation process, and adds further evidence for the number of 
singlet states contributing to the optical cycle of NV centres. Recently an infrared (IR) 
emission channel was observed, and attributed to a dipole-allowed transition between 
two singlet levels with an energy splitting of AE = 1.189 eV [231 [24] . In the same 




Figure 3. Decay of population from singlet states leads to recovery of fluorescence: 
(a) individual traces and (b) full data set of NV fluorescence with two excitation 
pulses and variable delay, (c) Fluorescence counts integrated over first 30 ns of second 
pulse as function of inter-pulse delay. The exponential increase is caused by decay out 
of singlet states, (d) Temperature dependence of singlet decay rate. The fit assumes 
a phonon-assisted decay process. 



publications it was shown that the IR emission follows the same time dependence as 
the visible transition, implying a short lifetime of the upper singlet state. Therefore 
the previously observed long-lived singlet state was attributed to the ground state of 
this IR transition. The temperature dependence of this singlet state lifetime has been 
reported for an ensemble of NV centres based on IR absorption measurements [24]. 
Here, we present data obtained on single NV centres. Since the oscillator strength of 
the IR transition is very weak, we use an indirect way to obtain this timescale, following 
the method used in [20]: We first apply a green pulse (1/xs) to reach a steady state 
population in the singlet states. This population manifests itself in reduced steady-state 
fluorescence with respect to the beginning of the pulse. After a variable delay we apply 
a second 1 /is pulse [figure E(a,b)]. During the delay between the pulses, population 
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stored in the singlet states decays back to the triplet ground state. The fluorescence at 
the beginning of the second pulse is proportional to the population in the triplet ground 
state and thus we can attribute its dependence on inter-pulse delay to the population 
decay out of the singlet state. In figure |3]^c) we integrate the photons emitted during 
the first 30 ns of the second pulse for each inter-pulse delay and fit these data to an 
exponential function. Figure EJ^d) summarizes the timescales we obtained in this way, 
for temperatures ranging from T = 13 K to T = 300 K for three different NV centres. 

We model the lifetime r of the singlet state as a combination of temperature- 
independent spontaneous decay rate Tq 1 and a rate accounting for stimulated emission 
of phonons of energy AE with an occupation given by Bose-Einstein statistics: r = 
To[l — exp(— AE/ksT)]. The fit yields a spontaneous emission lifetime of r = 371 ±6ns 
and a phonon energy of AE = 16.6 ± 0.9 meV, in reasonable agreement with ensemble 
data obtained in [21] . 

This result suggests that, apart from the 1.189 eV splitting, a third level AE below 
the IR transition's ground state is involved in the optical cycle of the NV centre. We can 
exclude that this third level is the triplet ground state, just 16.6 meV below the lowest 
singlet, as this would imply phonon-assisted spin relaxation on a sub- microsecond time 
scale at room temperature. This scenario clearly contradicts experimentally observed 
spin lifetimes on a millisecond timescale [5j[8]. Consequently the third level is likely to 
be another singlet state. The presence of three singlet states in between the triplet 3 A2 
and 3 E states was recently predicted by an ab-initio calculation of the excited states in 
the NV centre [25J, however, there a larger energy splitting between the lowest singlet 
states was obtained. 

6. Polarization probability 

Polarizing the electron spin by off-resonant optical excitation is a key technique for 
room-temperature spin manipulation of NV centres. Although this effect already has 
been identified to be caused by a spin-dependent ISC rate [20j [311 [32j [33], little is 
known about the relative contributions of spin-flip transitions between triplet states 
(fc 3 2, &4i), and ISC rates (fc 35 , fc 45 ,/c 51 , fc 52 ). We address this question by determining 
the polarization change due to a single excitation cycle. 

For that purpose we first initialize the NV spin by a 2 /is polarization pulse. After 
a waiting time of 1 /xs we excite the NV centre by a reference ps-pulse and measure 
the spin polarization by analyzing the relative contributions of the amplitudes in a bi- 
exponential fit to the fluorescence decay curve, as outlined in the previous section. This 
polarization corresponds to the steady-state value after CW excitation. To determine 
the change in polarization per excitation cycle we now apply a MW pulse to transfer the 
m s = population into the m s = ±1 states and then drive individual excitation cycles 
by applying 10 consecutive ps-pulses separated by 2 /is. For each excitation cycle, we 
again determine the spin polarization [figure HJ^c)]. All these bi-exponential fits use the 
same two time constants, obtained from a fit to the sum of all decay curves. From a 
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Figure 4. (a) We first polarize the spin by means of a 2 /is laser pulse at A = 532 nm. 
After a delay of 1 (is we apply a sequence of 11 pulses of 62 ps duration at A = 532 nm. 
For NV J, we apply a MW pulse 1 /is after the first ps-pulse to invert the spin state (for 
NV C we alternately run sequences with and without MW pulse). For each ps-pulse 
we measure the time-resolved emission. The experiment was performed at T = 300 K. 
(b) Fluorescence decay for consecutive ps-excitation pulses. The decay follows a bi- 
exponential function. A reference ps-pulse after a green 2 (is spin polarization pulse 
yields the initial optically induced spin polarization. A subsequent MW pulse transfers 
polarization to |2). (c) Change in polarization between consecutive pulses yields spin- 
flip probabilities Pi2(|l) — > |2)) and p2i(|2) — > |1)) per optical excitation cycle. 



power dependence measurement of the NV fluorescence rate we determine an excitation 
probablity a = 0.95 ± 0.05. 

The effect of a single excitation cycle on the spin polarization can be described 
by two counter-acting probabilities: a spin-flip from |1) to \2) (pi 2 ) and the opposite 
process (p2i)- Here, we consider only optically induced effects, i.e. the time between 
excitation pulses A t is assumed to be much shorter than the spin-lattice relaxation time 
(this assumption is substantiated by the constant polarization Pes for the 10 consecutive 
reference pulses in case of NV C [figure H^c)]). For state \m), population just before 
pulse n is denoted by P m , n , and just after the excitation pulse by P' m n . The steady state 
value P m , n=00 is abbreviated by P m . The pulse separation A t is much larger than the 
singlet decay time, such that P\. n + P2,n = 1 and therefore Pas,n = P\,n- 

Experimentally, we obtain spin polarization in the excited state. This differs from 
the ground state polarization due to a fraction e = ^23/^13 = ^14/^24 of spin non- 
conserving transitions. The populations for the (n+ l) th excitation pulse are then given 
by: 
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Pl, n+ 1 = ap 2 lP 2 ,n + (1 - aPl2)P\,n (1) 
P 2 ,n+1 = CtPYlP\,n + (1 - ap 2 l)P2,n (2) 

^3,n+i = a \T+ e - p 2,n+i + YjT^ Pl ' n+ y ^ 

Pin+l = « {^/W + Y^ P 2,n+l) • (4) 

The asymptotic value of the polarization is Pes = Vi\l{pvi + P21), and the excited 
state P ES , n = PU( P 3,n + K, n ) = [ P i,n (1 - e) + e] / (1 + e) « P 1>n closely follows the 
ground state spin polarization for small e. From [20] we can obtain an upper bound 
of e < 0.04. The change in polarization per pulse [figure BJc)] can be fitted by an 
exponential function Pss(n) = Pes + aexp(— n/c). From the steady-state polarization 
Pes an d the polarization rate c we can extract the spin-flip probabilities p 2 \ and pi 2 per 
optical cycle. These probabilities depend only weakly on e. Results are summarized in 
table m 

Table 1. Summary of parameters (T = 300 K), taking a = 0.95 ± 0.05 and 
e = 0.01 ± 0.01. p\2 and P21 are the total spin-flip probabilities per optical cycle. 
P35i P45 ; P5i an d P52 are the spin-dependent ISC probabilities for population in state 
|3), |4> and |5). 



II NV •' 


NV C 


Pl2 


0.078 ± 0.002 


0.079 ± 0.004 


P21 


0.315 ± 0.011 


0.372 ± 0.017 


Tigris) 


13.26 ± 0.03 


13.1 ± 0.1 


Ti A (ns) 


6.89 ± 0.06 


7.0 ± 0.2 


P35 


0.14 ± 0.02 


0.17 ± 0.03 


P45 


0.55 ± 0.01 


0.56 ± 0.02 


P5l/P52 


1.15 ± 0.05 


1.6 ± 0.4 



If we assume only rates indicated in figure HJb), use experimentally obtained 
lifetimes T 13 , T 14 and take e = 0.01 ± 0.01 and a = 0.95 ± 0.05, we can obtain the spin- 
dependent ISC probabilities of P45 = ^45/(^41 + ^42 + ^45) andp35 = ^35/(^31 + ^32 + ^35), 
and for the reverse process P51/P52 = ^51/^52; as summarized in table [TJ The low values 
of P51/P52 are in contrast to the current picture of the singlet relaxation process, where 
P52 = was used [20J. This gives another hint that an additional l E state needs to 
be considered for the relaxation process. The temperature dependence of p^i/p^ could 
give further insight into this topic. 

We note that this set of parameters can be used to model the spin- and power 
dependence of time-resolved NV centre emission using /xs excitation pulses at A = 532 nm 
within the framework of this five-level model, with reasonable agreement to experimental 
data [figure [5]. 
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Figure 5. Time- resolved excitation-power dependence of NV fluorescence excited by 
off- resonant 1 /is pulses (experimental data and model). The first pulse displays the 
fluorescence of NV centres that are optically polarized into m s = 0. The second pulse 
is obtained after applying a MW 7r-pulse, representing the fluorescence of a NV centre 
polarized in m s = ±1. The model shows the population in states |3) and |4). The 
calculations are based on parameters from NV J, assuming e = 0.01 and a — 1 (no fit 
is applied). Excitation rates are taken as integer multiples of 4 MHz. 

7. Summary 

We have experimentally determined the spin-dependent lifetime of the NV centre's 
excited state, whose difference is dominated by a spin-dependent ISC rate. Knowledge 
of these lifetimes allows us to determine the degree of spin polarization. In a second 
experiment, we identified the total lifetime of the singlet states, and, by analyzing its 
temperature dependence, the energy splitting of the long lived singlet transition. The 
measured energy of ~ 16meV indicates that at least three singlet states are involved 
in the optical cycle of the NV centre. Finally, we determined the spin-dependent ISC 
probabilities by analyzing the change of spin polarization induced by a single excitation 
cycle, without making assumptions on number and nature of the singlet states. The 
obtained rates are consistent with spin-dependent NV fluorescence dynamics based on 
a five- level model. 
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